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We studied the magnetization reversal rates of thermally assisted spin transfer torque switching 
in a synthetic free layer theoretically. By solving the Fokker-Planck equation, we obtained the 
analytical expression of the switching probability for both the weak and the strong coupling limit. 
We found that the thermal stability is proportional to Ao(l — I/I c ) 2 , not Ao(l — I /la) as argued 
by Koch et al. [Phys. Rev. Lett. 92, 088302 (2004)], where I and I c are the electric current 
and the critical current of spin transfer torque switching at absolute zero temperature, respectively. 
The difference in the exponent of (1 — I /Ic) leads to a significant underestimation of the thermal 
stability Ao. We also found that fast switching is achieved by choosing the appropriate direction of 
the applied field. 
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I. INTRODUCTION 

Spin transfer torque switching of the magnetization in 
ferromagnetic nanostructures has been extensively stud- 
ied both theoretically 1 " - — and experimentally 4- — because 
of its potential application to spin-clectronics devices 
such as magnetic random access memory. For device ap- 
plications, a thermal stability Ao = M H an V/ (2ksT) of 
more than 40 is required to guarantee retention time of 
longer than ten years, where M, H an , V, k-Q and T are 
the magnetization, the anisotropy field, the volume of the 
free layer, the Boltzmann constant, and the temperature, 
respectively. 

Recently, Hayakawa et a/.— showed that the 
anti-ferromagnetically coupled synthetic free layer, 
CoFeB(2.6nm)/Ru(0.8nm)/CoFeB(2.6nm), in a 
CoFeB(fixed layer)/MgO/CoFeB/Ru/CoFeB mag- 
netic tunnel junction shows a large thermal stability 
(Ao > 80) compared to a single free layer. On the other 
hand, Yakata et al%*^ showed that a ferromagnet- 
ically coupled CoFeB/Ru/CoFeB synthetic free layer 
shows a large thermal stability (Ao = 146 ± 29 for 
CoFeB(2nm)/Ru(1.5nm)/CoFeB(2nm) and 248 ± 60 for 
CoFeB(2nm)/Ru(1.5nm)/CoFeB(4nm)) compared to the 
single and the anti-ferromagnetically coupled synthetic 
free layer. These intriguing results spurred us to study 
a thermally assisted spin transfer torque switching in 
synthetic free layer. In contrast to the large number of 
experimental studies^—, few theoretical studies have 
been reported. Although the analytical expression 
of the switching rate of the thermally assisted spin 
transfer torque switching for the single free layer— ~—, 
P = 1 - exp[-/ t exp{-A (l - I/I c )(l - H appl /H an ) 2 }}, 
has been widely used to fit the experiments [see Eqs. 
(l)-(3) in Refs.— where I c is the critical current 
of the spin transfer torque switching at absolute zero 
temperature, it is not clear whether this single layer 
formula has validity when applied to a synthetic free 
layer. Thus, it is important to derive an analytical 



expression of the switching rate of the thermally assisted 
spin transfer torque switching for the synthetic free 
layer. 

In this paper, we studied the thermally assisted spin 
transfer torque switching rate for a synthetic free layer 
by solving the Fokker-Planck equation. The analyti- 
cal expressions of the switching rate were obtained for 
weak and strong coupling limits of the Fi and F2 lay- 
ers. One of the main findings was that the dependence 
of the thermal stability A on the current / is given by 
A oc Ao(l — I/I c ) 2 , not (1 — I/Ic), as argued by the pre- 
vious authors: — We emphasize that even for the single 
free layer A is proportional to (1 — I/I c ) 2 . The differ- 
ence in the exponent of the factor (1 — I/I c ) leads to a 
significant underestimation of the thermal stability Ao. 
We found that in the presence of the applied field -ff apP i, 
the switching times of the anti-ferromagnetically and the 
ferromagnetically coupled synthetic layers are different, 
and that fast switching is achieved by choosing an appro- 
priate direction of i? a ppi- 

This paper is organized as follows. In Sec. [TTJ we 
introduce the Fokker-Planck equation for the synthetic 
free layer and its steady state solution. We also introduce 
approximations to obtain the analytical expression of the 
switching probability. In Sees. IIIII and IIV1 we present 
the calculation of the switching probability in the limits 
of the weak and the strong coupling of the Fi and F 2 
layers. In Sec. [Vj we compare our results with those of 
other works. Section [VTl summarizes our findings. 



II. FOKKER-PLANCK EQUATION FOR A 
SYNTHETIC FREE LAYER 

Let us first derive the Fokker-Planck equation for the 
synthetic free layer. The system we consider is schemat- 
ically shown in Fig. [1] The two ferromagnetic layers, 
Fi and F2, consist of a synthetic free layer with the 
coupling energy —JSmi ■ 1112. Here m/j = M^/Mfc = 
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FIG. 1: The schematic view of the synthetic free layer con- 
sisting of the Fi and F2 layers, and p are the unit vec- 
tors along the directions of the magnetizations in the Ft, and 
the fixed layers, respectively, and d k is the thickness of the 
Ffc layer. i/ apP i, flan, and J represent the applied field, the 
anisotropy field, and the coupling between the Fi and F2 
layers, respectively. The flow of the electrons along the +x 
direction corresponds to the negative electric current I < 0. 



(sin 9 k cos ifk , sin 9 k sin tp k , cos 9k ) is the unit vector along 
the direction of the magnetization Mt of the layer. J 
and S are the coupling energy per unit area and the cross- 
sectional area of the system, respectively. It should be 
noted that J > and J < correspond to the ferromag- 
netically coupled and antiferromagnetically coupled syn- 
thetic free layers, respectively. Although we consider the 
ferromagnetically coupled system below, our formalism 
is applicable to the antiferromagnetically coupled system 
by changing the sign of the coupling constant J. We as- 
sume the uniaxial anisotropy along the z axis for both Fi 
and F2 layers, and the magnetizations mi and m 2 point 
to the positive z direction in the initial states. We also 
assume that the external field H app \ is applied along the 
z axis. Then, the total free energy F of the Fi and F2 
layers are given by 

F = - MiH^xVx cosf?! - -Mrifa^I/! cos 2 0j 

- M 2 H appl V 2 cos 2 - ^M 2 H an2 V 2 cos 2 9 2 

+ 27rA/ 2 V r i(sinf5i cos^i) 2 + 2nM 2 V 2 (sin9 2 cos ip 2 ) 2 

— JS [sin 9\ sin 9 2 cos(<pi — ip 2 ) + cos 9\ cos 9 2 ] , 

(1) 

where H ank , Vu — Sdk and dk are the uni-axial 
anisotropy field, the volume and the thickness of the 
Fk layer, respectively. The fifth and sixth terms in Eq. 
([T]) represent the magnetic energy due to the demag- 



netization field. We assume that \H. 



appl| 



to 



+e z to mi,ni2 = — e z . Following Brownii, we use the 
Fokker-Planck equation approach to calculate the switch- 
ing probability per unit time, where the Fokker-Planck 
equation is derived from the equations of the motion of 
the magnetizations. 

We assume that the dynamics of the magnetizations 
of the Fi and F2 layers are described by the Landau- 
Lifshitz-Gilbert (LLG) equations. In general, the spin 
transfer torque acting on mi arises from the spin currents 
injected from the fixed layer and the F2 layer. However, 
in a conventional synthetic free layer, the spacer layer 
between the Fi and F 2 layers consists of Ru, whose spin 
diffusion length is comparable to its thickness^; thus, 
the spin current injected from the F2 layer is negligible 16 . 
Then, the LLG equation of mi is given by 



dmi 



7imi x Hi + 7iajmi x (p x mi) 
dmi 



(2) 



71 mi x hi + aimi x 



di 



Similarly, the spin current injected from the Fi layer into 
the F2 layer is also negligible, and the LLG equation of 
m 2 is given by 

dm 2 dm 2 

- 7 2 m 2 x H 2 - 72m 2 x h 2 + a 2 m 2 x 



di 



di ' 



(3) 



where j k and a k are the gyromagnetic ratio and the 
Gilbert damping constant of the F^ layer, respectively. 
The magnetic field acting on the magnetization m^ 
is defined by = — (M k V k )~ l dF / dm k . h. k represents 
the random field on the ¥ k layer whose Cartesian com- 
ponents hki (i = x, y, z) satisfy 



(hki(t)h k/j (t')) 



2k B Ta k 
lkM k V k 



Skk'SijSit-t'), (4) 



guarantee at least two local minima of the free energy. 
When Hj k < -ffanfc, the states (mi,m 2 ) = (e z ,e 2 ), 
(e z ,—e z ), (— e z ,e z ), and (— e Zl — e z ) correspond to the 
energy minima, where Hj k = J/(M k dk). On the other 
hand, when H Jk > H &nkl the states (mi,m 2 ) = (e z ,e z ) 
and (— e z , — e z ) correspond to the energy minima. 

The purpose of this paper is to investigate the switch- 
ing rate of the magnetizations mi and m 2 from mi , m 2 = 



where (• • • ) means the ensemble average. Here we as- 
sume no correlation between the random fields acting on 
the Fi and F2 layers. The aj = Kr)l '/ \2eM\V\) term in 
Eq. ([2]) represents the spin transfer torque due to the 
injection of the spin current from the fixed layer. Here 
I is the electric current flowing along the x axis. The 
positive electric current corresponds to the electron flow 
along the — x direction, rj is the spin polarization of the 
electric current which characterizes the strength of the 
spin transfer torque. The explicit form of 77 depends on 
the theoretical model ^ 17 ' 18 , and, in general, depends on 
1 . However, for simplicity, we assume that ry is constant 
(the dependence of 77 on 9\ can be taken into account by 
replacing ajcosfli in Eq. ((6]) with J dcosf?ioj). p is the 
unit vector along the direction of the magnetization of 
the fixed layer. 

From the LLG equations @ and ©, we obtain the 
Fokker-Planck equation for the probability distribution 
of the directions of the magnetizations, W(m.i,m 2 ), 
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FIG. 2: The dependences of the effective potential & \v\, U2 
on 81 and 82 for (a) the weak coupling limit and (b) the strong 
coupling limit. The values of the parameters are written in 
Sees. IHII and HVl The white arrows indicate the most probable 
paths of the switchings mi, m.2 = +e z —¥ —e z . 



magnetization m^- from the yz plane increases the mag- 
netic energy due to the demagnetization field, it is rea- 
sonable to assume that the most probable reversal pro- 
cess is the magnetization reversal in the yz plane. In 
this limit, the demagnetization field plays no role on the 
calculation of the switching probability. By fixing the 
values of (pi and 922 to ir/2 or 3ir/2, the steady-state 
solution of the Fokker-Planck equation © is given by 
Wq oc exp[— ^/{ksT)], where the effective free energy 
is given by 



& = F - 
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which is given byll 
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Here we approximate that 1 + a k ~ 1 by assuming that 
afe <C l 19 . We also neglect the term proportional to 
aaj by assuming that \aj\ < |Hfc|, which is valid in the 
thermally assisted switching region. 

As shown by Brown-Li, the switching rate of the single 
ferromagnetic layer without spin transfer torque can be 
derived by using the steady-state solution of the Fokker- 
Planck equation and the continuity equation of the par- 
ticles of an ensemble [see Sec. 4. C in Ref. 14 ]. In the 
case of two ferromagnetic layers, as considered in this 
paper, the switching is described by the particle flow 
in (01,951,02,992) four-dimensional phase space, and, in 
general, it is very difficult to obtain an analytical ex- 
pression of the switching rate because the particle flow 
in the phase space is very complicated. To simplify the 
problem, we use the following two approximations. 

First, we assume that the magnetization rotates in the 
yz plane during the switching. Since the deviation of the 



The switching probability is calculated by using Wq oc 
exp[-J7(fc B T)]. 

The second approximation is that we consider the 
switching for the weak and strong coupling limits, where 
the weak (strong) coupling means that the magnitude of 
the coupling energy of the Fi and F2 layers, | — JSmi ■ 
ni2 1 , is much smaller (larger) than the uniaxial anisotropy 
ajMil^ sin 0i J W energy MkH an kVk/2. In other words, the weak (strong) 
coupling limit corresponds to i? an fc 3> ( < C)-ffjfc- The 
weak and strong coupling limit can be realized by chang- 
ing the thickness of the nonmagnetic layer g?n between 
Fi and F 2 layers because the magnitude of the coupling 
constant J strongly depends on c^n [see, for example, 
Fig. 2 in Ref4 or Fig. 1 (c) in Refj 1 ^], and varies from 
H j ~ 1 X 10 3 Oe to 1 Oe or less. Figures [2] (a) and (b) 
show the dependences of & on (0i, 02) for the weak and 
the strong coupling limit, respectively, where the white 
arrows indicate the most probable paths of the switching 
(the values of the parameters are written in Sees. IHII and 
IIVI with aj/a c i = 0.5). In the weak coupling limit, the 
magnetization reversal is divided into two steps: First 
mi reverses its direction from mi = +e z to 1112 = — e z 
by the thermally assisted spin transfer torque effect while 
the direction of rri2 is fixed to m2 = +e z , and second, 1112 
reverses its direction by the thermal effect and the cou- 
pling with the Fi layer while mi is fixed to mi = — e 2 . 
On the other hand, in the strong coupling limit, mi and 
ni2 reverse their directions simultaneously. 

By using the above two approximations, the calcula- 
tion of the switching rate is reduced to a one-dimensional 
problem. For the weak coupling limit, first we calculate 
the particle flow in 0i space, and second, we calculate 
the particle flow in 02 space. On the other hand, for the 
strong coupling limit, we calculate the particle flow along 
the direction of 0i =02- For such one dimensional prob- 
lems, the calculation method developed by Brown 1 ^ is ap- 
plicable to obtain the switching rate with some revisions. 
In Sees. IIIII and [TVl we show the switching probabilities 
for the weak and strong coupling limit, respectively. 

At the end of this section, we give a brief comment on 
the first approximation. The influence of the first approx- 
imation is that the critical current density estimated in 
our calculation, I c oc (i? a ppi + H an ), does not include the 
effect of the demagnetization field 4ttM [see Eqs. (|15p . 
([To]) . (|2"2"|) and (|2"3")) ] while the critical current density esti- 



sin 02 dif2 



(5) 
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mated by the LLG equation includes the demagnetization 
field, that is, J c oc (H app \ + H &n + 2it M) [see, for example, 
Eq. (14) in ReiM\. Since A%M > |ff ap pi|, #an, the criti- 
cal current density in our formula (10 5-6 A/cm 2 ) is much 
smaller than the experimental values (10 6 ~ 7 A/cm 2 )^— . 
One way to solve this discrepancy is as follows. Although 
we consider the in-plane magnetized system, it should 
be noted that our calculation is directly applicable to 
the perpendicularly magnetized system where the sys- 
tem has uni-axial symmetry and the switchings in the 
weak and strong coupling limits are described by only 9\ 
and 02- Suzuki et al. 20 showed that the effect of the de- 
magnetization field on the switching rate of the in-plane 
magnetized system can be taken into account by replac- 
ing iJan in the switching formula of the perpendicularly 
magnetized system by H an + 2nM . By applying this re- 
placement to our formula, our formula may be applicable 
to analyze the experiments quantitatively. The validity 
of this replacement requires the numerical calculation of 
the Fokker-Planck equation, and it is beyond the scope 
of this paper. 



III. WEAK COUPLING LIMIT 



In this section, we derive the switching rate of the mag- 
netizations for the weak coupling limit (H an k 3> Hjk) 
[see also the Appendix]. With this limit, the magneti- 
zation reversal is divided into two steps, as mentioned 
in Sec. |TTJ For convenience, we label the three re- 
gions around the potential minimum in the phase space, 
(#1,6*2) = (0, 0), (7T, 0), and (n,Tr), as regions 1, 2, and 
3, respectively. The first step (mi reverses from +e 2 to 
— e z ) corresponds to the transition of the particle from 
region 1 to region 2 while the second step (ni2 reverses 
from +e 2 to — e z ) corresponds to the transition from re- 
gion 2 to region 3. 

The switching rate from region 1 to re- 
gion 2 is obtained as follows 1 ^. In regions 
1 and 2, the distribution W(#i, 02 ) is given 
by W 1 exp[-{^(0 1 ,O) - J?(0, 0)}/(k B T)} and 
W 2 exp[-{^(0 1 ,O) - J?(7r,0)}/(fc B T)], respec- 
tively, where W x = W(0,0) and W 2 = W(ir,Q). 
The numbers of particles in region 1, n\, is ob- 
tained by integrating W(0i,O) over [0,# m i], where 
ml = cos _1 [-(iJ a p pl + Hji + aj/aij/H^x] gives the 
local maximum of the effective potential J^(#i,0). The 
explicit form of n\ is given by n\ = 2T4 7 ie^°'°)/( feBT ) ii, 
where factor 2 arises from the fact that we re- 
strict the particle flow in the yz plane; that is, 
ipi = 7r/2 or 3"7r/2 (in the anisotropic system con- 
sidered by Brown 1 ^, the numerical factor is 2tt, not 
2, as shown in Eq. (4.26) of Ref. 14 ). The integral 
h = J 9ml d6>isin(9iexp[-jF((9i,0)/(fcBT)] can be 



approximated 
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(7) 



The numbers of particle in region 2, ri2 = 
2W / 2e' 5r( - 7r:0 - ) / ( - feBT - ) /2, is obtained in a similar way by re- 
placing the factors &(0, 0) and <9 2 J^O, 0)/<96> 2 to J?(tt, 0) 
and d 2 ^ (tt , 0) / d0 2 , respectively. Next, we consider the 
particle flow from region 1 to region 2, I\^2- From the 
Fokker-Planck equation (JSJ, the particle flow along the 
0i-axis, Jg 1 , which satisfies li_>2 — 2sin#iJ6i 1 , is identi- 
fied as 



J 9i 



'M x Vi 



dF 



ajMiVi 



sin 6>i W + k B T 



dW 
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By multiplying e ^( g ^o)/(k B T) to 7^/(2 sin 6»i) = J 6l 
and integrating it over [0,7r], we find that [(712//2) — 
(n 1 //i)]/2 = -[M 1 V l /{2a lll k B T)}I 1 ^ 2 Irai, where the 
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The relation between the particle numbers in region 2 
and 3, 712 and n^, and the particle flow from region 2 
to region 3, i2->3j is obtained in a similar way. Then, 
by using the continuity equations of the particle flow, 
ni = -ii->2, n 2 = h^2 - h^3, and h 3 = 7 2 ^3, we 
find that the transitions of the magnetization directions 
among the three states, (0i,02) — (0, 0), (n, 0), (tt, it), are 
described by the following differential equations: 




-V12 V21 
^12 -(^21 + ^23) ^32 
V2Z -v Z2 




(10) 



The switching probability per unit time from the region 
i to the region j is given by = fij exp(— A^), where 
the attempt frequency and the thermal stability Ay- 
are, respectively, given by 
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Here J/J cl = a, 7 /a cl and J/I c2 = aj/a c2 . a cl (a c2 ) is the 
critical spin-transfer torque field to induce the magneti- 
zation reversal from region 1 (2) to region 2 (1) at zero 
temperature, and their explicit forms are given by 

a c i = -ax (H app i + Hjx + H anl ) , (15) 



a c2 = Oil (--ff a ppl - Hji + H anl ) , (16) 

respectively. Since |-ff a ppi + -Hji| is assumed to be smaller 
than i? an i, we find that a c i < and a c2 > 0. It should be 
noted that the description of the transition of the mag- 
netization by Eq. (fTO]) is valid for \a.j\ < \a c k\ because 
if \aj\ ^> |a c i | (|ctc2 1), the point mi = +e 2 (— e z ) would 
be unstable, and then we could not discuss the thermally 
assisted transition. We also note that the switching prob- 
abilities of m 2 , v 23 and ^32, are reduced to those obtained 
by Brownie by omitting Hj 2 where v 23 and 1/32 are in- 
dependent of the current /. 

When / is nearly Id, we find that ^12/^21 ~ 
exp[M 1 H anl V 1 /(2k B T)} > 1. Similarly when -iT apP l + 
Hj2 > 0, we find that ^23/^32 ~ exp[2M 2 (-i? a ppi + 
H J2 )V 2 /{2kBT)\ > 1. Within these limits, the ana- 
lytical solutions of Eq. (TTU|) with the initial conditions 
ni(0) = 1, 722(0) = 0, and 77 3 (0) = 0, are given by 
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Equation p^|) is the central result of this section: It com- 
pletely describes the magnetization switching of the syn- 
thetic free layer within the weak coupling limit. 

Figure |3] (a) shows a typical time evolution of 711 (t) , 
n 2 (t), and 773(7;) for a synthetic free layer with M = 995 
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FIG. 3: (a) Time evolution of n\(t), n,2 (t), and «3(t) with 
I /Id = 0.7 for the weak coupling limit, (b) The dependence 
of the switching rate 713(4) on the ratio I /Id for the weak cou- 
pling limit. For I /Id > 0.8, the switching time is saturated. 
The horizontal axis is the logarithmic scale. 



emu/c.c, H, 
7 = 1.732 x 
nm 
Fi= 



n = 50 Oe, ffappi = Oe, a = 0.007, 

10 7 Hz/Oe, d = 2 nm, S = tt x 70 x 160 

, and T = 300 K (for simplicity, we assume that 
F2 ^£ 1 io 1 2i > The current is taken tQ be = q 7 



The coupling constant is assumed to be J = 5.0 x 10~ 3 
erg/cm 2 , which corresponds to Hj = 25 Oe. From Eqs. 
([T7| , (fT8|) and ([T9"l) , one can easily see that the time evo- 
lution shown in Fig. [3] (a) is determined by two time 
scales, and v% 3 , which correspond to the switching 
rates of the Fi and F2 layers, respectively. Figure [3] (b) 
shows the dependence of the switching rate 773(7;) on the 
ratio I /I c i. For the currents |/| > 0.7|/ c i|, the switching 
times are on the same order (10 ms for our parameters), 
and for the large currents |/| > 0.8|J c i|, the switching 
times are saturated. This is because the current deter- 
mines the switching time of the Fi layer only, and for a 
large current, the total switching time of mi and m 2 is 
mainly determined by the switching time of m 2 , which 
is independent of the current. We can verify the satura- 
tion of the switching time from Eq. (|19p . where v\ 2 be- 
comes much larger than 1*23 as I approaches Id and then, 
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which is independent of the current i\ 



On the other hand, in the low current region \I\ <C \Id\> 
V12 becomes comparable or smaller than ^23, which leads 



6 



1.0 



5 0.5 





0.001 







/ - 


/0.8 


A 7 


/|/I C 1=0.6 



0.01 



0.1 1.0 
time [ms] 



10 



100 



FIG. 4: The dependence of the switching rate ri2(t) on the 
ratio I /Id for the strong coupling limit. The horizontal axis 
is the logarithmic scale. 



and is expressed as h\ = —h 2 — —v\ 2 n\ + ^21^2- The 
switching probability i/y = /y exp(— Ay) is given by 
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(20) 



to 713(f) ~ 1 — e - " 12 '. Then the switching time strongly 
depends on the current value because the switching time 
of mi becomes important to the total switching time. 
For example, the switching time for I /I c \ = 0.6 is longer 
than 100 ms, as shown in Fig. [3](b). 

The dependence of the switching time on the coupling 
constant J is as follow. By increasing the magnitude of 
Hj, the switching time rapidly decreases because of the 
fast reversal of the magentization of the F2 layer. For 
example, for Hj — 40 Oe with I/I c = 0.9, the switching 
time is on the order of 10 ~ 2 ms, which is three orders 
of magnitude faster than that for H j — 25 Oe. On the 
other hand, by decreasing the magnitude of Hj, only the 
magnetization of the Fi layer reverses its direction while 
the magnetization of the F2 layer remains m2 = +e z . 
For example, for Hj = 5 Oe with I/I c = 0.9, the switch- 
ing time of the Fi layer is on the order of 10 -2 ms while 
the switching rate of the F2 layer, n^, is approximately 
zero (n 3 ~ 10~ 9 ). Since the switching of the F 2 layer is 
induced by the coupling with the Fi layer, it is required 
to increase the magnitude of the coupling constant J for 
the fast switching by using a thin nonmagnetic spacer, al- 
though the increase of the coupling constant leads to the 
increase of the magnitude of the critical current density. 
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where I / I c i = aj/a c \ and I / I c2 = aj/a c2 , and the criti- 
cal spin-transfer torque fields in the strong coupling limit 
a ck are given by 



del = -on 



M 2 V 2 

H &pp \ + H an \ + (-ffappl + H an2 j 



a C 2 = OLi 



M x Vi 

M 2 V 2 
M-xVi 



(22) 



(— H app \ + H ari2 ) 



.(23) 

respectively. The analytical solutions of the transition 
equations, fi\ = —h 2 = —v\ 2 ri\ + ^21^2, with the initial 
conditions ni(0) = 1 and 712(0) = are given by 



rai(t) 



1/21 



V\2 



V\2 + V 2 \ V\ 2 + U21 



(24) 



IV. STRONG COUPLING LIMIT 

In this section, we derive the switching rate of the mag- 
netizations for the strong coupling limit {H ank -c H Jk ). 
For this limit, instead of (^1,^2) phase space, it is conve- 
nient to describe the particle flow in (*f?,i/j) phase space, 
where \& = (#1 + 2 )/2 and ip = 8± — 6 2 - Since mi and 
m 2 reverse their directions simultaneously, the reversal 
is described by the particle flow along the 'I'-axis with 
■0 = 0. For convenience, we label the two regions around 
the potential minimum in the phase space, ('J, ip) = (0, 0) 
and (71", 0), as regions 1 and 2, respectively. The conti- 
nuity equation of the particle in the regions 1 and 2 is 
obtained in a way similar to that described in Sec. IIIII 



„ 2 ( t ) = ^2_ e - (Kia + Wl )t. (25) 

V\2 + v 2\ V\2 + v 2\ 

When I is nearly I cl , 1^2/^21 ~ exp[(Mi/faniVi + 
M 2 H an2 V 2 )/(2k B T)} > 1. For this limit, Eqs. and 
(|25|) are reduced to 

m^-e-^', (26) 

n 2 (t) ~ l-e^ 12 *. (27) 

Equation l|27p is the central result of this section: It com- 
pletely describes the magnetization switching of the syn- 
thetic free layer within the strong coupling limit. 
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(a) (b) 




time [ms] time [ms] 

FIG. 5: (a), (b): The schematic views of the alignments of the 
magnetizations mi and m2 of (a) ferromagnetically (F) cou- 
pled and (b) anti-ferromagnetically (AF) coupled synthetic 
free layers after mi reverses its direction from the initial state 
(mi = +e z ) to mi = — e 2 . The directions of the applied field 
-Happi (assumed to be positive) and the coupling field H j2 are 
also denoted, (c), (d): The time evolution of 713(f) of the 
F-coupled (ng F '(i)) and the AF-coupled (n^ ¥ \t)) synthetic 
free layer with (c) H app \ — +5 Oe and (d) H app \ — —5 Oe. 
The value of the current is taken to be I/I c i = 0.7. 



For the strong coupling limit, the switching time 
strongly depends on the current I for all current region. 
Figure[4]shows the dependence of n 2 (t) on the ratio I/I c i, 
where the parameters used are same as those in Fig. [3] 
except J. The coupling constant J is assumed to be 
5.0 x 10 -2 erg/cm -2 , which corresponds to Hj = 250 
Oe. The orders of the switching times are 10 -2 ms for 
7/I c i = 0.8, 0.9, 1 ms for 7/I cl = 0.7, and more than 100 
ms for I / I c i < 0.6 in our parameter region, as shown in 
Fig. |U Such strong dependence of the switching time on 
the current arises from the thermal stability A12, which 
is proportional to (1 — I/Id) 2 , as shown in Eq. (|2"Tj) . 



V. RELATION TO OTHER WORKS 

In this section, we compare the results obtained in 
the previous sections to the other works^— . The top- 
ics discussed here are (1) the comparison of the switch- 
ing time of the ferromagnetically (F) and the anti- 
ferromagnetically (AF) coupled synthetic free layers, and 
(2) the comparison of the dependence of the thermal sta- 
bility to that obtained by Koch et alr^. 

First, we discuss the switching times of the F and 
the AF-coupled synthetic free layers. The difference 
in the switching times of these two kinds of synthetic 
free layer appears in the weak coupling limit with finite 
-ffappi- In this case, the switching time of the F-coupled 
synthetic free layer is characterized by Eqs. (|12p and 
(ITU) . For the AF-coupled synthetic free layer, the fac- 



tor +(-)(H app i - Hj 2 )/H a , n2 in Eq. $14$ is replaced by 
-(+)(i? ap pl + Hj 2 )/H an 2 while Eq. (JT5J) remains the 
same. This replacement is due to the fact that after mi 
reverses its direction from +e z to — e z , the sum of the ap- 
plied field -ff a ppi and the coupling field Hj 2 acting on m.2 
is -ff a ppi — Hj2 for the F-coupled synthetic free layer while 
it is -ffappi + Hj2 for the AF-coupled synthetic free layer, 
as schematically shown in Figs. [5] (a) and (b), and leads 
to the difference in the switching times of the F-coupled 
and the AF-coupled synthetic layers. 

The important point is that the fast switching is 
achieved by choosing the appropriate direction of H app \. 
Figures [5] (c) and (d) show the time evolutions of 713 (t) 
(Eq. ([TS])) for the F-coupled (n 3 F) ) and the AF-coupled 

(n 3 AF ^) synthetic free layers with (c) i7 app i = +5 Oe 
and (d) H app \ = —5 Oe. The current is taken to be 
I /Id = 0.7. The switching time of the AF (F) coupled 
synthetic free layer is faster compared to that of the F 
(AF) coupled synthetic free layer for H app \ > 0(< 0) be- 
cause both -ffappi and Hj2 assist the reversal of 1112. On 
the other hand, by changing the direction (sign) of -ff apP i, 
the switching time increases significantly because of the 
exponential dependence of the switching time (~ 1/f) 
on A cx [1 — (H app i±Hj)/H aD \ 2 . The difference between 

n 3 with -ffappi > and with H app \ < arises from 
the dependence of the switching time of mi on the direc- 
tion of -ffappi, and becomes negligible as I approaches Id 
because the total switching time is mainly determined by 
that of 1x12 within the limit of I /Id — >■ 1, as mentioned 
in Sec. IIIIl For the strong coupling limit, the switching 
times of the F-coupled and the AF-coupled synthetic free 
layers are the same because the coupling energy is con- 
stant during the switching in this limit, and the coupling 
field Hj plays no role on the switching. 

Second, we discuss the dependence of the thermal sta- 
bility A on the current /. As shown in Eqs. (fT2"j) and 
(|2"Tj) . our calculations show that A oc A (l — I / I c ) 2 ■ It 
should be noted that our formula is applicable to the sin- 
gle free layer by omitting the coupling of the Fi and the 
F2 layers, and thus, even for the single free layer we find 
that A oc Ao(l — I/I c ) 2 . Recently, a similar result was 
obtained by Suzuki et aZi^S and Butler et al£2- for the 
perpendicularly magnetized single free layer. However, 
the formula of the switching rate with A oc Aq(1 — I / I c ) 
first obtained by Koch et al*^ has been widely used to 
fit the experiments^—. 

The important point is that the difference of the ex- 
ponent of (1 — I/Ic) leads to a significant underesti- 
mation of Ao Let us consider the fit of the exper- 
imental results of the switching rate with the formula 
P = 1 — cxp{— f t exp[— Ao(l — I/I c ) n ]}, where for sim- 
plicity we assume that the attempt frequency is con- 
stant fo- When I/I c = 0.5, the thermal stability Ao es- 
timated by our formula (n = 2) is two times larger than 
that estimated by the conventional formula (n = 1). 

The difference between the exponent of (1 — J/7 C ) in 
our calculation and that in the theory of Koch et alr^ 
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arises from the steady-state solution of the Fokker-Planck 
equation of the free layer magnetization: 



dW 
dt 



1 8 



MV sin 9 89 
'8F 



ajMV 



sin ( 



W - 



, 8W] 



(28) 



Koch et at argued that the steady-state solution of 
of Eq. d2HD is W K och oc exp{-F[l + aj/(aH)]/(k B T)}, 
where H = |H| is the absolute value of the magnetic field 
acting on the free layer magnetization. However, when H 
depends on m, WkocIi is not a steady state solution of Eq. 
(|28p . In general, H depends on m because of the presence 
of the uni-axial anisotropy field H an = H an m z e z , which 
guarantees two local minima of the free energy F. Thus, 
in the calculation of the switching rate, we should use 
W oc exp{-[l - {a j MV cos 9) / (aF) / (k B T)}} , which is 
the steady state solution of Eq. (|2"5)l as shown in Sec. lU 
instead of WkocIi • The difference between Wo and WkocIi 
leads to that of the exponent of (1 — I / J c ) in A. 



VI. CONCLUSIONS 

In conclusion, we studied the magnetization switching 
of the synthetic free layer theoretically by solving the 
Fokker-Planck equation. We obtained the analytical ex- 
pression of the switching rate for the weak and the strong 
coupling limits, given by Eqs. (TT9|) and (j27|) . We found 
that the switching time within the weak coupling limit 
becomes saturated as the current / approaches the crit- 
ical current I c \. We compared the switching time of the 
ferromagnetically and the anti-ferromagnetically coupled 
synthetic free layers with a finite applied field, and find 
that fast switching is achieved by choosing the appro- 
priate direction of the applied field. We also found that 
the dependence of the thermal stability on the current is 
A cx A (l - I/Ic) 2 , not A oc A (l - I / I c ) as argued by 
previous authors^, which leads to a significant underes- 
timation of Aq. 
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APPENDIX A: DETAILS OF THE 
CALCULATION IN SEC. []Tj] 

In this appendix, we show the details of the deriva- 
tion of Eq. [see also Sec. 4 C in Ref^]. First, 
let us consider the switching from region 1 to region 2. 



The number of the particle in region 1 is obtained by 
integrating 2W X exp[-{J*"(0i, 0) - J?(0, 0)}/(fc B T)] over 
< m i < 0mi ; that is, 



ii = 2W 1 e^ ( -°' 0)/( - kBT) / d0isin0iexp 



^(0i,O) 



It should be noted that the exponential term in the in- 
tegral rapidly decreases by changing 0i from to m i- 
Then, we replace J?(0i,O) by its Taylor series about 
0i = 0, keep the terms up to the second order of 0i, 
and replace the upper limit of the integral by oo. The 
first term of Taylor series, 8^(0, 0)/90i, is zero because 
0i = corresponds to the local minimum of sin0i 
is approximated to 0\. Then, we arrive Eq. ([7]). The 
numer of the particle in region 2, n 2 , is obtained in a 
similar way; that is, n 2 = 2W 2 e^ ty ' K ' " > / tykBT " > I 2 , where I 2 
is given by 



k B T 



<9 2 J^(tt,O)/<902 



The particle flow from region 1 to region 2, I±. 
isfies [see Eq. ((SJ] 



8W 
~df 1 



I gg. 



W = - 
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MxVi 



ctijik B T ) 2 sin 0i 



(30) 
►2) sat- 

(31) 



According to Browni^, we assume that 1\^, 2 is indepen- 
dent of 0i. By multiplying e ^(e u o)/{k B T) to Eq _ ^ and 

integrating it over [0, tt], the left hand side of Eq. (f3~Tj) is 
reduced to 

o 90i 

= W(ir, o)e^fi)/(k B T) _ w ^ 0)e ^(o,o)/(fc B T) (32) 
_ 1 / n 2 ni 

where we use the definitions of 1\ and I 2l i.e., n\ — 
2ff ie ^(W)/(teT) /l and n2 = 2W 2 e^ { * fi)/{kBT) h. On 
the other hand, by using Taylor series of J^(0i, 0) about 
0i = 0mi, the right hand side of Eq. (f3"Tj) is approximated 
to -[MiVi/(2ai7ifc B T)]Ji^2-fn l i; where I ml is given by 
Eq. ©. Thus, we obtain 



ri2 
h 



ni 
h 



M x Vi 



2ai7i/c B T 



h- 



(33) 



Similarly, the number of the particles in regions 2 and 
3, n 2 and 713, and the particle flow from the region 2 to 
region 3, /2->3, satisfy 



"3 n 2 



I' 3 



I2 



M 2 V 2 



<2^3-<m2, 



2a 2 ~j 2 k B T 

where I 2 , I3, and J m 2 are, respectively, given by 

k B T 



J 1 = e -^(7T,0)/(teT) 



8 2 ^(Tr,0)/89l 



(34) 



(35) 
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(36) 



^23 



V M 2 V 2 J I 2 I m2 ' 



(40) 



Im2 = 



2irk B T 



3 ^(ir,e m2 )/(fcBT) 



d^(7r,6 m2 )/d9 2 2 



sin y m2 



(37) 



where 6» m2 = cos ^-(flappl - H J2 )/H an2 }. 

By using Eqs. (j33]), (f3"4"j) and the continuity equations 
of the particle flows, Ui = — Ii-^j "2 = ^1^2 — -^2^3, 
and n3 = 72->3, we obtain Eq. (ITU1) . The switching 
probabilities per unit time, i/y , are given by 



^12 



^21 



«i7ifcBJ 1 
MiVx 



1 



hi 



ai7ifc B r 



l-'ml 



f 



(38) 



(39) 



^32 



/ «272 



V M 2 V 2 



(41) 



The explicit forms of Vij are obtained by using ^(61, 62) 
and its derivative, and given by Eqs. (|lll) - (fl4")) . 

By applying the negative current and magnetic field, 
which induce the magnetization reversal from mi , m 2 — 
+e z to m 1 ,m 2 = — e z , v^ 2 and ^23 become much larger 
than v 2 \ and ^32, respectively. Then, Eq. (|10l) can be 
approximated to h\ = —v\ 2 ri\, h 2 = v\ 2 ri\ — v 2 ^n 2 and 
fi3 = v 2 j,n 2 . Then, the solutions of ni, n 2l and 713 with 
the initial conditions, ni(0) = 1, n 2 (0) — 0, and n 3 (0) = 
0, are given by Eqs. (|T7)) - (fTO)) . 
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